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Abstract. Properties of solutions of generic hyperbolic systems with multiple 
characteristics with diagonalizable principal part are investigated. Solutions are 
represented as a Picard series with terms in the form of iterated Fourier integral 
operators. It is shown that this series is an asymptotic expansion with respect to 
smoothness under quite general geometric conditions. Propagation of singularities 
and sharp regularity properties of solutions are obtained. Results are applied to 
establish regularity estimates for scalar weakly hyperbolic equations with involutive 
characteristics. They are also applied to derive the first and second terms of spectral 
asymptotics for the corresponding elliptic systems. 

MSG 2000 classification: 35S30, 35L45, 35L30, 35C20, 58J40 

1. Introduction 

Let X be a smooth manifold without boundary of dimension n > 3. Let P be an 
elliptic self-adjoint pseudo-differential operator of order one acting on half-densities 
on m-dimensional crossections of vector bundles on X. We consider the following 
Cauchy problem for u = u{t, x) 



XI) 



iu' -Pu = 0, {t, x) G [0, T] X X, 
u\t=o = Uq. 



It is well known that if equation (jl.ll) is strictly hyperbolic, the system can be diag- 
onalized and its solution can be given as a sum of Fourier integral operators applied 
to Cauchy data (e.g. [Hj). An important question that has been studied over many 
years is what happens when P has multiple characteristics. 

Since we will be mostly interested in local properties of solutions, we may already 
assume that P acts on functions, and can think of it as an m x m matrix of pseudo- 
differential operators of order one and we think of uq as of an m-vector. 

Let A{x,^) denote the principal symbol of P. If A is a diagonal matrix, properties 
of system (jl.ip have been studied by many authors. For example, in [16^ and _15j 
Kumano-go and coauthors used the calculus of Fourier integral operators with multi- 
phases to show that the Cauchy problem is well-posed in L^, Sobolev spaces 
if*, and to study its propagation of singularities. Systems with symmetric principal 
part A have been extensively studied as well (e.g. [14j, |il2|, etc.) In a generic 
situation, such systems have double characteristics, and their normal forms have 
been found by Braam and Duistermaat j2] . Recently, Colin de Verdiere jlj used these 
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representations to derive some asymptotic properties of such systems. Polarization 
properties of similar systems were studied by Dencker in pj. 

More elaborate analysis of system (jl.lll becomes possible if we assume that the 
principal symbol matrix A{x,C,) is smoothly microlocally diagonalizable with smooth 
eigenvalues aj{x,^) and smooth eigenspaces. Then, as it was pointed out by Rozen- 
blum in [21], there exists a finite dimensional cover X of X such that A lifted to 
X can be globally digonalized provided X is compact. In this situation Rozenblum 
showed that the Picard series for this problem gives an expansion with respect to 
smoothness in the case of non-involutive characteristics. In other words, one assumes 
that if aj{x,^) = ak{x,^) for j ^ k, then the Poisson bracket {a^, afc}(x, ^) 7^ 0. This 
means that at all points of multiplicity, bicharacteristics intersect characteristics sur- 
faces transversally. However, this condition is non-generic even for diagonal systems. 
For example, it is clear that if one of characteristics has a fold, there may be a point 
where this transversality condition fails, and it is not possible to remove it by small 
perturbations. 

One purpose of this paper is to present results removing the transversal intersection 
condition. We will allow characteristics to be involutive of finite type and some 
characteristics to be identically equal. Operators satisfying our Condition C below 
will be also generic in the class of smoothly microlocally diagonalizable systems. 
Below we will explain that the microlocal diagonalizability condition is quite natural 
when considering weakly hyperbolic scalar equations with Levi conditions (Examples 
1 and 2). This is also the case for Maxwell equations (e.g. [1 ). 

We will investigate regularity properties of system (jl.ip in this generic setting. 
Even in the case when the system is strictly hyperbolic, properties of solutions 
have been studied for many years, since already this case has several important 
applications for nonlinear equations and harmonic analysis. Regularity properties of 
non-degenerate Fourier integral operators have been established by Seeger, Sogge and 
Stein in [21] . They showed that a Fourier integral operator T of order zero satisfying 
local graph condition, is locally bounded from (-L^)comp to Lf^^ for 1 < p < 00 and 
a = (n — l)|l/p— 1/2|. As a consequence they showed that if system is strictly 
hyperbolic, there is a loss of a derivatives in L^*, i.e. Uq G implies u{t, ■) G . 
Moreover, if at least one of characteristic roots aj is elliptic, the loss of a derivatives 
is sharp. If none of a^-'s is elliptic, this result can be improved ([22])- 

Our Theorem 12 . 21 will establish a similar property for systems fll . Ij) with multiplic- 
ities. Moreover, this will imply estimates for scalar weakly hyperbolic equations 
with involutive characteristics. It is known that in general weakly hyperbolic cases 
one often loses regularity even in L^. However, in the case of involutive characteristics 
the equation can be diagonalized and in Theorem 12.31 we will give U' estimates for 
such equations. This will, on one hand, extend the U' result of Seeger, Sogge and 
Stein to weakly hyperbolic equations and systems with multiplicities, while on the 
other hand establishing U estimates for systems considered by Kumano-go, Rozen- 
blum, and others. The result will be general for scalar weakly hyperbolic equations 
satisfying Levi conditions with characteristics satisfying Condition C below. 

Note that if a scalar operator strictly hyperbolic and we write it in the form 
we can diagonalize P together with lower order terms (e.g. [T3]) and spht it into m 



2 



scalar equations, for which many things are known. However, in the case of multiple 
characteristics this is impossible, so a more elaborate analysis is needed. 
Now we will formulate our main assumption. Let us define operator 

HaJ = {ajJ}J = l,...,m, 

where 



n . 

L — 1 V 



dg_d£_ _ dg df 

d^k dxk dxk d^k 



k=l 

is the usual Poisson bracket. Our assumption is that at points of multiplicity aj = ak 
of non-identical characteristics aj and ak, bicharacteristics of aj intersect level sets 
{ttk = 1} with finite order, i.e. H^^ak 7^ for some M, at points where aj = ak- 
in other words, we allow involutive characteristics of finite type, and formulate our 
main condition. 



Condition C: 

Suppose that there exists M E N such that if for some j and k, aj and are not 
identically the same, then 

aj(a;,0 = afc(x,0, {x,OeT*X^ 

(1-2) .JU 

Haj^kix, = {aj, {aj, ■ ■ ■ Wj, ctfc }} • • •} {x, 0, 

for some A < M. While the function M = M{x,C,) is locally bounded, it is allowed 
to grow at infinity. 

We note here that the transversality assumption of Rozenblum |2I] requires (|1.2|) 
to hold with M = 1. Strictly hyperbolic case is also covered by this condition (in 
which case we set M = 0). The case of aj and defining glancing hypersurfaces (as 
in Melrose ^7|) corresponds to M = 2. 

In Section 2 we will give several examples of characteristics satisfying condition C, 
in particular those arising from weakly hyperbolic scalar equations with involutive 
characteristics. Such equations and propagation of their singularities have been ana- 
lyzed in [3], [HI, [HI; [U]; [in|; ^tc. We will also estabhsh estimates in and other 
spaces for the weakly hyperbolic equations or systems satisfying condition C. 

It is interesting to note that conditions similar to Condition C appeared in the 
study of subelliptic operators (e.g. Hormander [TUi Chapter 27]). For instance, in 
the case of 2 x 2 systems P with characteristics Oi and 02, we can consider operators 
Q with principal symbol q = ai + ia2- Then microfocal subellipticity of Q implies 
Condition C, with some M dependent on the loss of regularity for which, therefore, 
implies the Weyl formula for P f Theorem 12. 5|) . regularity estimates for and all 
other results of this paper. 

Now we will give an informal explanation of the strategy of our analysis. First, let 
us follow [21] to show that microfocal diagonalizability implies a local one on some 
cover X of X with finitely many leaves. For this argument we assume that X is 
compact. Since all the analysis of this paper will be local, if X is not compact, we 
can always assume that the amplitude of P{x, D) is compactly supported. 
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Let an elliptic pseudo-differential operator P{x, D) of order one act on sections of an 
m-dimensional Hermitian vector bundle E. Let L'^{E) be the space of sections of half- 
densities on E and let P be self-adjoint on L'^{E). Let E' be the lifting of E to T*X. 
Then for each (x, ^) G T*X the principal symbol A[x,^) of P{x,D) is a Hermitian 
isomorphism of sections of E'. Without loss of generality we can assume that A{x,^) is 
positive definite. Indeed, if it has both positive and negative eigenvalues, it is possible 
to globally block-diagonalize A{x, D) with some suitably chosen pseudo-differential 
operators, to reduce it to a direct sum of positive and negative definite operators. 
Then each of these operators can be analyzed independently. 

Let us assume that the principal symbol A{x,C,) is microlocally diagonalizable. 
This means that microlocally in A C T*X such that E'\\ = Ax C™^, principal 
symbol A{x,C,)\a has m smooth non-negative eigenvalues aj{x,C,) and one dimensional 
eigenspaces Vj{x,$,), and such diagonalizations are compatible in intersecting cones. 
In this situation Lemma 15.11 insures that there is a cover X of X with finitely many 
leaves such that the principal symbol of the lifting of P{x, D) to T*X can be globally 
diagonalized modulo lower order terms. Note that since dimensions of X and X are 
the same and because of formula (j5.1|) all our results on X will imply corresponding 
results on X. Therefore, we may assume that the principal symbol A of operator P 
may be smoothly diagonalized over compact subsets of X, that is 

P = A + B,A = diag{Ai, . . . , Am}, 

where Aj G \E'^ are scalar pseudo-differential operators with principal symbols aj(x, ^) 
(eigenvalues of A). Here a/s may be identically equal to each other or may intersect 
with any finite order, according to our Condition C. Here S is an m x m matrix of 
pseudo-differential operators or order zero. We may also assume that Bjj = for 
1 < j < m, if we add these terms to the diagonal of A. For the moment we will also 
assume that none of Oj's are identical. Otherwise, if some of Oj's being identically 
equal to each other locally at some points, the construction is slightly different, but 
all the results remain valid. This will be carried out in detail in Section 3 in the proof 
of Theorem 12.21 Substitution U = e~*^*\^ leads to the equation 



;i.3) 



V = Z{t)V, 
V\t=o = /, 



with Z{t) = —ie'^^^Be ^'^^ Writing the Picard series for problem ()1.3|) . we obtain the 
expansion 

(1-4) V{t) = / + £z(ti)dti + /o7J^Z(ti)Z(t2)«i + --- 

A general term of this series is 



I 



i-1 



(1.5) Qi= / ■■■/ e'^nH,,j,e'^^'^^''-'''^---dti...dti. 

Jo Jo Jo 

It is easy to see that | IQ;] |l2_^^2 < C/U and that series p.4j) converges in and 
in H''. Using the notion of a multi-phase for Fourier integral operators, Kumano-go 
et al. ([ini, [in]) studied propagation of singularities of Qi. Instead of introducing 
multi-phases for Qi, we will analyze operators Qi in more detail and will show its 
smoothing properties in Sobolev spaces under Condition C. Here, contrary to the 
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transversal non-involutive case of Rozenblum (when M = 1), we do not have good 
control on the singular supports of integral kernels of operators Qi, so more elaborate 
geometric analysis is required. 

In fact, our Theorem 12.11 asserts that Qi{t) maps to some Sobolev space Hp^^\ 
and that — > oo as / — oo. Then we will show that this allows to treat the Picard 
expansion as a series with finitely many terms, with many conclusions, such as V 
estimates for solutions and spectral asymptotics of P. 

In particular, estimates will follow from a general principle which we will prove 
in Theorem 13. II for equation u' — Z(t)u = f,u{0) = Uq. In Corollaries 13 . 21 and 13 . 31 we 
will show regularity of solutions of this equation for pseudo-differential operators of 
order zero Z{t) G or for Fourier integral operators of negative orders Z{t) G /~^. 
In general, there may be problems with this construction for Z{t) G /°, but there we 
can use the structural properties of Z{t) in (jl.3|) . 

Everywhere in this paper = '^^^{X) will denote the space of classical pseudo- 
differential operators of order fi of type (1, 0). The space of Fourier integral operators 
of order fi with amplitudes of type (1, 0) will be denoted by I^. All Fourier integral 
operators in this paper will be non-degenerate, which means that its canonical relation 
satisfies the local graph condition, i.e. it is a graph of a symplectic diffeomorphism 
from T*X to itself. Constants C may be different in different formulas throughout 
this paper. We will use the following notation for norms and spaces. By we will 
denote the Sobolev space of functions / such that (/ — A)"/^/ G L^. For a function 
/ we will denote its L^-norm by and its Sobolev if^ = norm by ||/||s- If T 

is an operator, by \\T\\s we will denote its operator norm from to H^. 

Results of this paper can be established also in the case of operators P dependent 
on t. This will be the subject of a separate paper since it will also involve the analysis 
of non-smooth coefficients. 

We would like to thank Ari Laptev for drawing our attention to Rozenblum's 
paper pi], Grigori Rozenblum for several valuable references and Yuri Safarov for 
discussions. The research was supported by EPSRC grant GR/R67583/01. 

2. Main results 

Now we will formulate our results concerning terms of the Picard series (jl.4j) and 
solutions to systems (jl.lj) and (jl.3p . We will also give a Weyl formula for P. 

Our first main result will be Theorem 12.11 on the smoothing properties of terms 
Qi. The other two important results will be Theorem 12.21 on the L^-regularity of 
solutions to system ()1.1|) and Theorem 12.51 on the spectral asymptotics for elliptic 
operator P satisfying Condition C. To obtain L^-estimates, we use Theorem 13. H 
which we regard as a general principle behind regularity estimates for general Cauchy 
problems based on several natural properties of the right hand side operators Z{t). 
We will illustrate its use in several situations in Corollaries 1-3. Theorem 12.41 is a 
statement on the propagation of singularities of operators Qi or, more generally, of 
Fourier integral operators in which the frequency integration is performed over a cone 
rather than over the whole space. Note that these results will hold for microlocally 
smoothly block-diagonalizable operators with any (finite) geometry of characteristics, 
i.e. characteristics satisfying our Condition C. 
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The following theorem establishes a smoothing property of operators Qi under 
Condition C. 

Theorem 2.1. Let condition C be satisfied, that is assume that there is M such that 
for any (x, ^) G T*X and any 1 < j, k < m, j ^ k, with aj and Ofc not identically the 
same and aj{x,^) = ak{x,^), we have 

X 

{aj, {aj, ■ ■ ■ {aj, Ofc }} . . .} (x, 7^ 

for some number X < M. Then for sufficiently large I operator Qi in is bounded 
from Ll^^p to Hf^^, where 

^ _ ((-3n/2-2)(3 [//2] -l-n) + {[1/2] -n- 1){1/{2M) - n - I)) I 



(3[//2] -2n-2 + //(2M)) 6M + 2 



Note that the exact order can be improved. However, it is most important that 
it increases to infinity as / ^ oo. This theorem implies, in particular, that the series 
()1.4p . i.e. the series 

V{t) = I + Qi{t)+Q2{t) + --- 

is a series with respect to smoothness. This fact allows one to refine the study of 
propagation of singularities and regularity properties of solutions to systems 
and (fO|l . 

It turns out that smoothing properties of Qi in Theorem 12. II can be used to establish 
local properties of solutions to systems with multiplicities (ll.lj) . For strictly 
hyperbolic equations such estimates have been established by Seeger, Sogge and Stein 
in [21], and some optimal estimates were given in [221 • O^'^ next result concerns 
regularity of solutions to Cauchy problem In this theorem we also allow the 

lower order term B to depend on t. 

Theorem 2.2. Let 1 < p < oo and a = {n - l)\l/p - 1/2|. Let P = P{t,x,D^) be 
an m X m matrix of elliptic classical pseudo-differential operators of order one. Let 

Pit, X, D^) = A{x, D,) + B{t, X, a,), 

where A is a symmetric matrix of pseudo- differential operators of order one and 
B is a matrix of operators of order zero. Assume that the matrix A is smoothly 
(microlocally) diagonalizable, with smooth eigenspaces and real eigenvalues aj{x,C,), 
satisfying condition C. Then for any compactly supported f G L^H L"^, the solution 
u = u{t, x) of the Cauchy problem 

du 

(2.1) ^_-P{^t,x,D,)u = Q, uiO) = f, 

satisfies u{t, ■) G Lf^^ for all < t < T. Moreover, there is a constant C > such 
that 

sup ||M(t,-)||L. <CT||/|Ug. 
0<t<T 
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We note that it is sufficient to only assume that A is microlocally diagonalizable. 
Then by Lemma 15.11 we can first have the statement of Theorem 12.21 on the cover 
X. Then, using formula ()5.ip and the fact that dimX = dimX, we get the same 
conclusion on X. 

As a consequence, if the Cauchy data Uq is compactly supported, we obtain local 
estimates in other spaces as well: 

• Uq E L^s+a implies u{t^ ■) G L^, s G M. 

• Uo G Lip(s + {n — l)/2) implies u{t, ■) G Lip(s). 

• Let 1 < p < g < 2. Then uq G -^^^_i/g+„/p_(„_i)/2 implies u{t,-) G L^. Dual 
result holds for 2 < p < g < oo. 

The proof of Theorem 12.21 is based on the other main result Theorem 13.11 which we 
will discuss in the next section. Estimates in other spaces follow by standard methods 
of harmonic analysis ( j2Sl ) • 

Now we will give some examples where condition C holds while the tranversal- 
ity condition (M = 1) fails. This is, for example, the case when pairs aj,ak define 
glancing hypersurfaces or when we consider Maxwell systems with variable coeffi- 
cients. Below we will concentrate on systems arising from scalar weakly hyperbolic 
equations with Levi conditions. 

Example 1. In scalar equations with Levi conditions studied by Chazarain 
Mizohata-Ohya [THl, Zeman j^, one assumed that {a^, a^} = Cjk{aj —ak). It is clear 
that in this situation aj{x,^) = ak{x,^) implies {aj,ak}{x,^) = 0. However, in a 
general case when Cjk{x,C,) is non-constant, condition C is satisfied generically. 

Example 2. Let L be a scalar operator with involutive characteristics. More pre- 
cisely, let us denote dj = Dt + \j(t, x, D^) and let 

(2.2) L = di---d^ + Y^ j^dj, ■ ■ ■ dj^ + c, 

k<m 

where b{t,x, Dx),c{t,x, D^) G \I'° are pseudo-differential operators of order zero for 
all t G [0, T]. We will assume that symbols of all operators are infinitely differentiable 
with respect to t in the topology of symbols of the corresponding order. Let us assume 
that operator L has involutive characteristics, i.e. that 

[dj, dk] = djdk - dkdj = ajkdj + Pjkdk + jjk, 

where ajk, Pjk,1jk G \E'° are pseudo-differential operators of order zero. Then it was 
shown by Morimoto in that the Cauchy problem for the equation Lu = / is 
diagonalizable (with 1 + Y2^=i ''^^J^- components). Even in the simplest case of 
characteristics not depending on t, we have 

{Xj, Xk} = ajk{x,^){Xj - Xk) + jjk, 

similar to Example 1. 

Propagation of singularities of systems with vanishing Poission brackets has been 
also studied in these situations. For example, Iwasaki and Morimoto jTHj studied 
propagation of singularities of 3 x 3 systems, where the second Poisson bracket vanish. 
Also, Ichinose JJJ studied 2x2 systems with vanishing second Poisson brackets. 
Theorem 12.21 implies a precise statement on estimates. 
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Theorem 2.3. Let 1 < p < oo, a = (n — l)\l/p — l/2\ and s G M. Let L be as in 
(j2.2p and suppose that principal symbols aj{x,^) of Xj satisfy condition C and do not 
depend on t. Let u be a solution to the Cauchy problem 



(2.3) 



Lu = 0, 

^m(0, x) = gjix), <j <m - 1, 

and let Cauchy data gj G L^_j^^ be compactly supported. Then u(t, ■) G {Lp)ioc for 
all t G [0, T] and 

m—l 

(2-4) sup \\d\{t,-)U<CY,\\g,\\,.^ ^, 

te[o,T] ^+^+™ 

where = dj^ . . . dj,^, k < m—l, and {ji, . . . being permutations of some elements 
o/{l,...,m}. 

Note that in the strictly hyperbohc weh as in some very special cases of 

operator L in ()2.2|) (e.g. when all b and c are zero), following the method described 
by Treves in and estimates for Fourier integral operators, it is possible to obtain 
the estimate for the Sobolev norm ||M||rp in the left hand side of (12.41) . 

" " s+m— 1 ^ ^ 

Let us now give a final example of L^ estimates for second order equations, which 
we will prove in the next section. 

Example 3. Let us consider the second order equation 

u" + b{x, Dx)u' + c(x, D.ji:)u = 0, 

where 6 G \E'\ c G Let us denote (x) = (1 + x^)^/^. Introducing v = ^ ^^^^ ^ , 
the matrix form of this equation is given by 

{D^) 

Let bi and C2 be principal symbols of b and c. The equation is hyperbolic if bf > 4c2, 
with multiple roots at 6f = 4c2. Assume that 6f — 4c2 = /i^, with fi & being a 
symbol of order one. Then characteristics 01,02 satisfy {01,02} = |{&i,/i}, which 
may vanish. 

Let the Cauchy data be m(0) = /o,m'(0) = fi. Let a = (n — l)|l/p — 1/2| and 
1 < p < 00. If /i is elliptic, it is known that if fj G L^_j, then u(t,-) G L^. Our result 
of Theorem 12.21 will imply that the same is true for any fi E S^. 

Note that in this example we may requre only microlocal diagonalization with the 
same conclusion. 

Let us now discuss the propagation of singularities for operators Qi. This result is 
essentially a reformulation of Rozenblum's result in the case of finite geometry under 
Condition C. It is clear (also from multi-phase analysis) that singularities propagate 
along broken Hamiltonian flows. Let 

J = {ji, ■ ■ -Ji+i}, 1 < jfe < m, jk ^ jk+i- 

Let ^j{t,x,^) be the corresponding broken Hamiltonian flow. It means that points 
follow bicharacteristics of Oj^ until meeting the characteristic of o^^ , and then continue 
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along the bicharacteristic of a^^, etc. Note that singularities may accumulate if wave 
front sets for different broken trajectories project to the same point of X. 
We can write 



Qi= [ I{t)dt, 

J A 



where i = {ti, . . . ,ti), A = {0 < < ti^i < . . . < ti < t} is a symplex in M} and 
= Z{ti) o . . . o Z{ti). It is possible to treat it as a standard Fourier integral 
operator with the change of variables t = Cl^l""*^- Let be a cone in = R"+'. Let 

Iu{x)= [ [ e'^^'''y^^^aix,y,9)uiy)dyde 
Jk Jy 

be a Fourier integral operator with integration over the cone K with respect to 6. Let 
Kj be or a face of K. Let (pj{x, y, 6j) = (p\K,,Oj e Kj. Let C T*X x T*X be a 
Lagrangian manifold with boundary: 

For G C T*Y, let Aj{G) = {z e T*X : 3C G G : (z, C) G AJ. Then we have the 
following statement on the propagation of singularities. 

Theorem 2.4. Let u G V'iJ). Then WF{Iu) C \^jkj{^ F{u)). 

The proof is standard and follows Hormander [9^ . 

From this, we can deduce first and second terms of the spectral asymptotic of 
operator P. Let us call T a period of symbol A{x^ ^) if there exists J such that 
ji = ji+i, and the trajectory of $j is closed: $j(T, x, ^) = (x, ^). Then we have the 
following extension of well-known results of Hormander [S], Duistermaat-Guillemin 
|7), Safarov-Vassiliev [2S], and Rozenblum pTj . 

Theorem 2.5. Assume that X is compact and assume that Condition C is satisfied. 
Let D he the set of (x, ^) G T*X such that there exist T and J such that ^j{T, x, ^) = 
(x, ^). Assume that the measure of D is zero. Then for the spectrum of P the following 
Weyl formula holds: 

iV(A) = : A, < A} = c„A" + <A"^i + o(A"-i), 

where Xj are eigenvalues of P. 

Proof of the L^ estimates will be based on Theorem 13.11 which we regard as an 
independent result on regularity of solutions of partial differential equations. Theorem 
15.31 concerns the measure of the set where a function is small given some information 
on the multiplicity of its roots. It will play a crucial role in the proof of the smoothing 
property in Theorem 12.11 



3. Regularity of solutions 

In this section we will present a principle governing solutions of first order systems. 
Let Z{t) G C{C^{X),V'{X)), t G [0,T], be a time dependent family of operators. 
Let Wo, Wi and W be hnear subspaces of V'{X) such that Wq, Wi ^ W. We will 
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make different choices of these spaces in the future, dependent on the structure of 
operators Z{t). Let us consider the Cauchy problem for u = u{t, x): 

. . j u'-Z{t)u = r, r{t) G Wq, 

^'^ > \ u{0) G Wi. 

One is often interested in the following question. If the right hand side and Cauchy 
data satisfy r(t) G Wq and m(0, •) G Wi, when do fixed time solutions u(t, ■) of (|3.ip 
belong to W7 In general, some loss of regularity is possible in problem (j3.ip even if 
Z{t) are very good. So we will think of Wo being the smallest, Wi an intermediate, 
and W the largest among these spaces. The following theorem says that if operators 
Z{t) have some structure, and solutions of Cauchy problem ()3.1|) with zero Cauchy 
data are in W, so will be solutions with Cauchy data from some sufficiently large 
space Wi. 

Theorem 3.1. Let Wo,Wi ^ W be linear subspaces ofV'{X). 
Let Z{t) G C{C^{X),V'{X)), t G [0,T]. Assume that 

(i) (Boundedness) Z extends to an operator in L°°{[0,T], C{L'^{X), L'^{X))) , and 
Z{t) extend to continuous linear operators from Wi to W , for all t G [0,T]. 

(ii) (Calculus) Z{ti) o ■ ■ ■ o Z(ti) : Wi W are continuous for all I and for all 
t^,...,ti G [0,T]. 

(iii) (Smoothing) There exists I such that 

z{t) [ r ... r ' z{ti) 0...0 z{ti)dti ...dti 

Jo Jo Jo 
is continuous from Wi to Wq, for all t G [0, T]. 

(iv) (Zero Cauchy data) Solutions v = v{t,x) of the Cauchy problem 

v' - Z{t)v = r, r{t) G Wo, 
v{0) = 0, 



(3.2) 



satisfy v{t, ■) eW forte [0, T]. 
Then the solution u = u{t,x) of the Cauchy problem 



(3.3) 



u' - Z{t)u = r, r{t) G Wq, 
u{0) G W,, 



satisfies u(t, ■) e W for all t G [0, T]. 

Moreover, if Wo,Wi,W are normed spaces and if solutions v{t,-) to ()3.2p in (iv) 
satisfy \\v(t, < ^^11^(^)1114^0 /^'^ '^^^ ^ ^ [O;^]; 'then also 

\\u{t,-)\\w<C{\\umw, + \\rmwo), 

for allte [0,T]. 

Conditions (i) and (ii) ensure that operators Z{t) have some structure. Indeed, if 
Wi C is different from W , (ii) does not follow from (i). In our typical applications, 
Z{t) will be time dependent pseudo-differential or Fourier integral operators, and 
compositions in (ii) are essentially of the form of a single operator Z{t). Condition 
(iii) is natural from the point of view of harmonic analysis, since integration with 
respect to a parameter often brings additional regularity. Condition (iv) ensures 
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that solutions with zero Cauchy data and regular right hand side are also sufficiently 
regular. 

Proof. Let U{t) be an operator solving the Cauchy problem 

j U'- Z{t)U = R{t), Rit) e CiW,, Wo), 
\ f/(0) = /. 

Let Uo(t) be some partial solution to the problem 

r U[,-Zit)Uoit) = R{t), 
1 Uo{0) = 0. 

Then the solution U of problem (j3.4|) satisfies 

(3.5) U{t) =Uo{t) + I + Z{ti)dti+ / Z{U)Z{t2)dt2dti + ... 

Jo Jo Jo 

The convergence of this series can be understood in L^. Indeed, because of assumption 
(i), the term of this series with k integrals can be estimated by sup^ 1 1^(^) Wl^^l'^/^^- 
From this it also follows that U{t) is a solution of ()3.4j) in L^. Let us now define 

SN{t) =1+ Z{ti)dti + Z{ti)Z{t2)dt2dt, + ... 

Jo Jo Jo 

+ / .../ Z{ti)Z{t2)...Z{tN)dtN ...dt2dti. 
Jo Jo Jo 

Let V{t) = U{t) — SN(t), it is equal to Uo(t) plus the remainder of the series ()3.5p . 
Then we have 

V'{t) - Z{t)V{t) = 
{U-SNy{t)-Z{t){U-SN)it) = 
iU'-ZU)it)-iS'j,-ZS^)it) = 

pt Pt2 ptN-l 

(3.6) R{t)-Z{t) / .../ Z{t2)...Z{tj^)dt2...dtN. 

Jo Jo Jo 

Choosing N = I, from assumption (iii) of the theorem the second term is continuous 
from Wi to Wq. Since also R(t) G C(Wi, Wq), it follows that the right hand side is a 
continuous linear operator from Wi to Wq. 

Let w = u{0) G Wi be the Cauchy data for ()3.3|) . If we denote by p(t) the value of 
the operator in the last line of ()3.6p at w, we will have p(t) G Wq. The value of V^(0) 
is 

1/(0) = U{0) - Sn{0) = 0. 

It follows now that V{t)w solves Cauchy problem ()3.2j) . so it belongs to W by as- 
sumption (iv). Since SN(t) is continuous from Wi to W by assumption (ii), and 
V{t)w = U{t)w - SN{t)w is in W, be obtain u{t, ■) = U{t)w G W. 

Moreover, suppose that we also have the estimate ||f (t, ■)\\w < C\\p{t)\\w^^ in (iv). 
Then we also have 

\Ht,-)\\w < \\V{t)w\\w+\\SN{t)w\\w < CMt)\\wo+C\\w\\w, < C{\\w\\w, + \\RitMwo)- 

□ 
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Later we will need this in the case of Z{t) being Fourier integral operators of 
order zero. However, let us point out several applications to other cases of pseudo- 
differential and Fourier integral operators. In these cases we will make different 
appropriate choices of spaces Wq, Wi, W. Moreover, in Corollaries . 21 and . HI we will 
assume that the corresponding non-homogeneous Cauchy problems with zero Cauchy 
data have unique solutions. This is a natural assumption if Z{t) behave sufficiently 
well with respect to t since we are working in subspaces of L^. 

Corollary 3.2. Let 1 < p < oo. Let Z{t) E t E [0,T], be a family of pseudo- 
differential operators of order zero with amplitudes compactly supported in x, y, uni- 
formly m t. Suppose that Z E L°°([0, T], C{L^, L^)) nL°°([0, T],C{Lp, 1^)). Then the 
solution u = u{t, x) of the Cauchy problem 

u' - Z{t)u = r, r{t) eLp, t> 0, 
u{0) E LP, 



(3.7) 



satisfies u(t, ■),u'(t, ■) E L^, for all t E (0,T]. Moreover, we have an estimate 

sup \\u{t,-)\\Lv < C\\u{0)\\lp. 

te[o,T] 

If amplitudes of Z(t) are not compactly supported with respect to x,y, we have 
a similar local statement for compactly supported Cauchy data. Note also that 
pseudo-differential operators of order zero are locally bounded in L^ and L^, for all 
1 < p < cx). Conditions Z E L°°([0, T], £(L2, L^)) and Z E L°°{[0,T], C{Lp, Lp)) 
simply mean that we have some control on their norms, i.e. there exist a constant C 
such that 

\\zit)\\L2_,L2 < c, <c,ytE [o,T]. 

Proof. Let us choose W = Wi = Wq = Lp^^^. Let us check conditions of Theorem 
13.11 Properties (i) and (ii) follow from regularity properties of pseudo-differential 
operators of order zero and our assumptions. Property (iii) also holds because Z{t) 
are locally bounded in Lp. Property (iv) is a consequence of Duhamel's principle and 
is similar to the one in Corollarv 13.31 Picard series is convergent in Lp provided that 
operators norms ||Z(t)||/^p^iP are uniformly bounded for t E [0,T]. Norm estimate 
follows from this as well. 

We can see that u'(t,-) E W from u' = Z{t)u + r and from the continuity of Z{t) 
in W. □ 

We will now apply Theorem 13. II in the case of Z{t) being Fourier integral operators. 
While our case ()1.3p corresponds to Z{t) being operators of order zero, the crucial 
smoothing property (iii) will follow from the fact that operators Z{t) have a special 
structure. For general Fourier integral operators Z{t) without structure, we need 
to assume that they are of negative orders. This is for example the case when the 
zero order term B in Theorem 12.21 is actually a pseudo-differential operator of some 
negative order. 

Corollary 3.3. Let I < p < oo, e > and a = {n - l)\l/p - 1/2|. Let Z{t) E 
t E [0,T], be a family of non- degenerate Fourier integral operator of order —e with 
amplitudes compactly supported in x,y, uniformly in t. Suppose that operators Z[t) 
can be composed and that Z E L°°{[0,T], C{H'^ , H'^)) , for some s > {pn — 2n)/2p 
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when p > 2 and s = when p < 2. Then the solution u = u{t,x) of the Cauchy 
problem 

u' - Z{t)u = r, r{t) e H', t> 0, 
u{0) e {LD n L\ 



(3. 



satisfies u{t, ■) G L^, for all t G (0,T]. Moreover, 

sup ||M(t,-)IUp<C||«(0)|Up. 

*e[o,T] 

Note that operators Z{t) are locally bounded in H"^, so assumption 
Z G L^{[0,T],C{H',H')) simply means that \\Z{t)\\Hs^H' < C for all t G [0,T]. 

Proof Let = n L^,^^, W^i = {Ll),,,^^, and P^o = ^^^^ C 1^. Let us check 
conditions of Theorem Kill Condition (i) follows from the fact that non-degenerate 
Fourier integral operators of order are bounded from {L'^^comp to Lf^^. Condition 
(ii) follows from the calculus of non-degenerate Fourier integral operators, since we 
assumed that compositions of Z{t) are again non-degenerate Fourier integral oper- 
ators. Smoothing condition (iii) for large / follows again from the calculus, since 
operators Z{t) are of order — e. 

Finally, let us show that solutions of v' — Z{t)v = r(t), r(t) G H^^^^, with zero 
Cauchy data f (0) = 0, satisfy v(t, ■) E L'^ H L^. In fact, we will show that v(t,-) G 

From the uniqueness of the solution of this problem it follows that we can use 
Duhamel's principle to write 

(3.9) v{t,x)= [ E{t,s)r{s,x)ds, 



where E{t, s) is the propagator of 

{dt-Z{t))E{t, s) =0, 



E{t, s)\t=s = L 

Picard series for this problem gives the asymptotic expansion of E{t,s), in partic- 
ular implying that E{t, s) is bounded in and in H'^ provided operator norms 
\\Z{t)\\H'>^H'> are uniformly bounded for t G [0,T]. From ()3.9|1 it follows that v(t, ■) G 
H^. Moreover, since r(s, ■) G H'', we also get an estimate 

\\uit, ■)\\h- < C sup ||r(r, ■)\\h^, 
^e[o,T] 

implying the estimate in Corollarv 13.31 □ 

Proof of Theorem \2.'/A As we have already mentioned, by Lemma f5. II we can assume 
that characteristics of A are correctly defined on T*X. Since A is diagonalizable, we 
can write 

(3.10) P{t,x,D^) = Qaj{x,D^) + B{t,x,D.,), 

Bit,x,D,) = iB,k{t,x,D,)),<,,k<m, B,keC^i[0,T],^°). 

Some of a/s may be identically equal to each other. We can renumber a/s into r 
groups (possibly of size one) of equal characteristics. These are the eigenvalues of the 
matrix A{x, ^) counted with multiplicity. Thus, we have 1 = ki < k2 < ■ ■ ■ < K = 
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n + 1, and a^- = . . . = ^ ctfc, for k < ki or k > fcj+i. This means that we 

have a group of the same roots ai, . . . , a^j-i, etc., while roots from different groups 
are not identically the same. Therefore, this is a decomposition of the first order 
principal part into a block-diagonal form with the same roots in each block, with 
possible equality of roots in different blocks at some points. So we can write 

(3.11) P(t,x,D^)= diag {di, ... ,ar) + B(t,x,D^), 

where di = diag (a^^, . . . , a^.^-^-i) are diagonal matrices with equal roots at the di- 
agonal. Let us set 

Ai = di + {Bf,^)k 

SO that 

p = A + B = diag (ii, . . . , A,.) + B. 
Note that in the last equality we can assume B^,j = for ki < fi,^ < fcj+i — 1 if we add 
these components to the corresponding components of A. Let U{t) = exp{—iAt)V{t). 
This is well defined in view of, for example, |26, VIII]. Then 

V = Z{t)V ^ -ie'^'Be-'^'V, V{0) = I. 

Now we will apply Theorem 13.11 with Z{t) = —ie^^^Be""^^^. Let us choose W = 
Llomp n LP, Wi = {LPjcomp, and Wo = i/^'omp with s > {pn - 2n)/2p for p > 2 and 
s = for 1 < p < 2. Conditions (i) and (ii) follow from the calculus and regular- 
ity properties of non-degenerate Fourier integral operators of order zero. Smoothing 
condition (iii) follows from Theorem 12.11 For condition (iv) we can use Duhamel's 
principle similar to the proof of Corollarv 13.31 Thus, Theorem 13. II implies that prop- 
agator V{t) is continuous from L^, to L^. Operator U(t) is a composition of V(t) 
with a non-degenerate Fourier integral operator exp{—iAt), so U{t) is given as a sum 
of a smoothing series obtained by the multiplication of Picard series for V{t) with 
exp{—iAt). Using the calculus of Fourier integral operators in each term of the series 
and its smoothing property we can repeat the argument of Theorem 13.11 in this case 
to see that u{t, ■) G with an estimate for its norm. □ 

Note that if B in Theorem 12.21 is a pseudo-differential operator of negative order, 
B G ^E^'^, for some fi < 0, the proof is simpler because we do not have to use Theorem 
12.11 to prove condition (iii) of Theorem 13.11 Instead, we can use directly Corollary 
13.31 to obtain the smoothing property (iii). 

Proof of Theorem \2.^ Let L be as in ()2.2|) and let u be the solution of 

(3.12) Lu = /, Diu{0, x) =9j{x), 0<j <m- 1. 
Let 

U = {u, diu, d2U, . . . , did2U, d2diu, . . . , d-^u, . . .)^, 
where J = {ji, . . . , jfc} is a permutation of some elements of {1, . . . , m}, \J\ = k < 
m — 1. Vector U has 1 + X^^^'Til/j! components. Here we can write d"^ = + 
J2j=oCj{t,x,D^)D{, where c/(t) e We set \J\ = k. It was shown by Morimoto 

in ^] that U solves the system 

(3.13) DtU + AU + BU = F, f/(0, x) = G{x), 
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where F = (0, 0, /,...,/) and G = {go, ...,g\j\ + T}j=o cjgj, . . .), ^ is a diagonal 
matrix with Aj's at the diagonal and 5 is a matrix of pseudo-differential operators 
of order zero. Matrix B has some operators in the last row, zeros, and —1 above the 
diagonal. If Xj satisfy Condition C, Theorem 12 . 21 implies \\U{t, ■)\\lp < C*||G'||j^p. Since 
Cjit) G we get < CY^"^=Q WQjWi^^ ^ which implies the estimate of 

the Theorem. □ 



4. Estimates for Picard series 

In this Section we will prove Theorem 12. li on the smoothing properties of terms Qi 
of the Picard series ()1.4j) . 

Let Aj G j = 1, . . . ,r, be elliptic pseudo-differential operators of order one. 
Let aj{x,^) denote their principle symbols. We can assume that there are no idential 
symbols among these a/s. Let 

H(t) = e*^w*ie*^M(*2-<i) . . . f^-i^3l+l^l ^ 
where 1 < jfc < r, ^ jk+i, k = 1, . . . ,1 + 1, and i = (ti...ti). Let us define 

Q= / ... / B{t}H{t}dti...dh, 
Jo Jo Jo 

where i?(t) G \E'° is a pseudo-differential operator of order zero smoothly dependent 
on t. Such operators Q appear in the Picard series (|1.4|) . (|1.5|) . In this section we 
will give a detailed description of operator Q in order to prove that it is a smoothing 
operator when / is sufficiently large. First of all let us note that if is a Fourier integral 
operator and due to the theorem on compositions of Fourier integral operators the 
canonical relation A* C T*X x T*X of H(t) is given by 

A* = {(x,p,y,e): (x,p) = **(l/,0}, 

where \E^* = <|)*.i o • ■ ■ <|)*'"*'"^ o $7*^ and is the Hamiltonian flow defined by a,-. 
It can be easily checked that H is a solution operator for the system of equations 

dH 

(4.1) — = iTk{h,...,tk)H, k = l,...,l, 

Otk 

where G ^E'^ is a pseudo-differential operator of order one. In view of Egorov's 
theorem its principle symbol is equal to 

(4.2) T°(ti, = (a,, - a,,^J o o • ■ • o (x, 0, (^,0 e T*X 

for all k = 1, . . . ,1. Let us construct a phase function (p{t,x,y,C,) which defines the 
operator H{t) for small i. We will look for it in the form ip{t, x,y,C,) = V'(t, x,C,)—y-C,- 
It follows from ()4.1|) and ()4.2j) that ip satisfies a system of Hamilton- Jacobi equations 



(4.3) || = T°(t,a;,|^), ^(O,a;,0 = x-e 

In [21] it was checked that Frobeneus conditions for system ()4.3|) are satisifed. Solving 
this system we obtain a non-degenerate phase function. This phase function defines 
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a Lagrangian manifold A , so that we have 

(4.4) ix,^) = ¥{y,0 = ix\y,0,p\y,0), y = ^- 

Now we are going to investigate the smoothing properties of operator Q. We can 
write Q as 

Qu{x)= [ [ [ e'^^''^'y'%{t,x,y,Ou{y)dyd^di, 

where if satisfies ()4.4j) and b is an amphtude of order zero, which we may assume to 
be compactly supported with respect to x and y. Here A is a symplex {0 < ti < 
ti-i < . . . < ti < t}. 

Let X ^ be a cut-off function such that xi"^) = 1 fo^' < 1 xi"^) = for 
|r| > 2. Operator Q can be decomposed as Q = Ri + R2 + R3, where 

R,u{x)= ! [ [ e'''^''^'y'^^fi,{e,i,x,y,Ou{y)dyd^di, J = 1,2,3, 

J A Jr" Jr" 



where 



f^2ie,i,x,y,0 = Xi^-'\^\)(l-xi^''\^m-')]bii,x,y,0, 



dC \ dt 
MeJ,x,y,0 = X{e''\^\)xie~'\^\\^nbit,x,y,0- 

Let us first consider operator Ri. On the support of fii we have the estimate 

dip 



> e. 



Therefore, there exists an operator L(-^) of order 1, such that Le**^ = e*"^, with coef- 
ficients estimated by Ce^^ on the support of the amplitude of fii. When integrating 
by parts with L there may appear an additional factor of when differentiating 
X- Integrating by parts p times with this operator L we obtain an operator with an 
amplitude of order —p and coefficients that can be estimated by e~^^. From Lemma 
15.21 with q = —p we obtain the following estimate 

(4.5) ||i?i||p_„_i<Ce-3p+"+i. 

The same procedure with integrating by parts with respect to ^ can not be applied 
to i?2- But here there is a possibility to integrate by parts with respect to i. Indeed, 
there exists an operator M(^), such that Me'"^ = e^'^, with coefficients not greater 
than e~^|^|~^ on the support of /i2. Integrating by parts with M we obtain an operator 
with an amplitude of order —1 with coefficients that can be estimated by e~^, where 
another may appear from differentiating x- The boundary integrals have the same 
form as Q but they have amplitudes of order —1 and depend on not less than / — 2 
time variables. The reason for possibly losing two variables is that after restriction 
to the boundary, say t2 = h, it may happen that and the same roots. It 

follows that we can apply such procedure [1/2] times. As a result we obtain operators 
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of order with coefficients not greater than Ce ^['/^l. Then in view of Lemma 

Ofor [1/2] > n + 1 we have 

(4.6) ||i?2||[z/2]-„-i<C6-3['/']+"+\ 

Let us now consider the last integral R^. It is not possible to apply procedures with 
integrating by parts as above either with respect to ^ or with respect to t. But in this 
case it turns out that the support of amplitude /X3 is small. The singular support of 
the integral kernel of Q may be very irregular in this case, so a more delicate analysis 
is necessary to show the smoothing properties of -R3. First we will show that 

(4.7) |i;°(t,xUo,/(i/,0)l <c^e|e|,j = i,...,/, 

on the support of /is. We notice that /is differs from only if 
(4.S) |||.., ,|,..|,. 

It follows from (jOI) and (gS)) that 

(4.9) |7;°(t,a;,|^)|<Ce|e|,j = l,...,/. 

Because of homogeneity of Tj* with respect to ^ we also have the following trivial 
estimates 

Consequently, we obtain 

dx 



\d^T^{t,x,^)\<C\^\, \d^T^{t,x,^)\<C, j = l,...,l. 



|i;°(t,a;,/)-T,°(t,x*(y,0,/(2/,0)l 



<Cm\x- x\y, 1 + 1 1^ - p\y, 01), J = 1, . . . , /. 



Since ||(t, a;*(?/, 0, Z/, = 0, we get 



^{t, X, y, = ^{t, X, y, - ^{t, x\y, ^),y, = -^^(^^ Oix - x\y, 0), 

q2 — 

for some points x*. Since | 7^ for small t, we obtain 



I dxd^ 

(4.10) \x-x\y,0\<C 



^it,x,y,i) 



From the properties of Lagrangian manifold A* in ()4.4|) we also obtain 

(4.11) |g(t,a;,0-/(y,OI = l/(^,0-/(y,OI < C|^(t,a;,y,Olle|, V(x,y,0. 

Finally, taking into consideration ()4.8p - ()4.11|) we obtain ()4.7p on the support of /is. 

Estimate ()4.7|) shows us that the support of the amphtude /is of operator i?s is 
contained in 

Si = {(t, X, : a: G X, 0, (t, y, ||) G S}, 
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where 

S = {{i, y, ^) : \T;(t, x\y, /(y, ^)) | < Ce, j = 1, . . . , I}. 

That allows us to apply estimate ()5.4j) of Lemma 15.21 with g = and 5 = 1 to the 
operator R3 to get 

(4.12) ||i?3||-3n/2-2 < Ce-"-imeas(H), 

where meas is the natural measure on [0,T]' x S*X. To use this inequality we need 
to estimate the measure meas(H). The set S maps to 

E2 = {{i,y,0-- \T;;ii,y,0\<Ce, \p\y,0\ = I, j = I, ■ ■ ■ ,1} 

under the Hamiltonian flow ^* which preserve measure. The measure of S2 may be 
estimated by the measure (in [0,T]' x X x R*^) of 

S3 = {(t,y,0: \T;;{i,y,0\<Ce, c<\^\<C, j = l,. ..,/}. 

Let us now introduce sets 

(ti, . . . , y, e) = {t, : \T!^{t, y,0\< Ce, c < \^\ < C}, j = 1, . . . , I. 

To estimate measures of these sets we will use Theorem 15.31 For this, in the notation 
of Theorem 15. 3^ we set w = (y,^) and consider the function 

t ^ fit, w) = T°(ti, . . . , tj-i,t, tj+i, ...,tuy,i), 

where t takes the position of tj in i. Condition on Poisson brackets shows us that 
function / may have zeros in t of order not greater then M. It follows now from 
Theorem 15.31 that 

(4.13) maxmeas{t,- G [0,T] : T^{i,y,0 < Ce} < Ce^/^M^ 

where meas is just the Lebesgue measure, y varies over a compact set and c < |^| < C. 
Now we can use a simple observation that if we have two functions /i(ti) and f2{ti, ^2) 
such that 



meas{ti : |/i(ti)| < e} < Ce°, maxmeas{t2 : 1/2(^1,^2)! < e} < Ce' 

ti 



then 



meas{(ti,t2) : 1/1 (ti)| < 6,1/2(^1,^2)1 <e} = / dhdt2 

\Mtl)\<e,\f2iti,t2)\<€ 

dt2 ] dti 

|/l(tl)|<e \J\f2{tut2)\<e 

Applying this argument / times to the estimate ()4.13j) . we get 

meas(S3) < C vol (supp 6) max 11 max meas(Sj(t, ?/, e)), 
which implies 

meas(S) < Ce'/'*'^. 
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Here we used that the support of the amphtude b{t,x,y,^) is compact with respect 
to X and y, and that T is smalL Finally, combining this with estimate (j4.12|) . we get 

(4.14) ||i?3||-3n/2-2<C6'/^'"-"-\ 

An application of Lemma EIEl to ()4.5|1 . ()4.()j) and ()4.14j) yields the requred estimate 
for Q in Theorem 12.11 

5. Various auxiliary results 

First we describe the lifting of the problem to insure that characteristic roots of 
the principal symbol of P are globally uniquely defined. The proof follows the paper 
of Rozenblum j21j. We give it here for the completeness and since we will need this 
construction to determine the orders for estimates in Theorem 12.21 

Lemma 5.1. Let X be a smooth compact manifold without boundary of dimension 
n>3. Then there exists a cover X of X with finitely many leaves, such that on the 
lifting E of E to T*X branches of eigenvalues aj{x,^) and eigenspaces Vj{x,^) of the 
principal symbol A{x,^) are smooth and globally well defined. The space L'^{E) has 
a decomposition into a direct sum of m spaces such that the matrix representation 
of P{x,C,) with respect to this decomposition consists of pseudo-differential operators, 
and its principal symbol is a diagonal matrix with aj{x,C,) at the diagonal. 

Proof. Let us fix zq = (xo,^o) ^ S*X. For each path in S*X beginning at Zq we look 
at the continuation of the diagonalization along this path. At each point z G S*X we 
obtain up to permutation several possible collections aj{z), Vj{z), such that homotopic 
paths from zq to z correspond to the same collection. Since for n > 3 homotopy of 
paths in S*X and their projections to X is equivalent, we obtain a homomorphism 
TTiiX) to the permutation group of order m. A construction in the homotopy theory 
( |20j ) gives a cover p : X X with finitely many leaves, such that the lifting of this 
homomorphism to vri(X) is trivial. So, each closed path in S*X takes eigenvalues 
and, therefore, also eigenspaces, to themselves. This means that eigenvalues and 
eigenspaces have global smooth branches on S*X and hence also on T*X. 

Let p^{x, ^) be a family of orthogonal projectors on Vj{x, By the standard Gram- 
Schmidt process, we can add lower order terms to p^{x,C,) to obtain symbols Pj{x, ^), 
for which pj{x, D)pk{x, D) = 6jkPj{x,D) and ^Pj(x,-D) = 1. For u G L^{E), let 
Uj = pj{x,D)u. Since P{x, D)pj{x, D) = pj{x, D)P{x, D)pj{x, D) modulo lower 
order terms, P = YliPj^Pk yields the desired diagonalization. □ 

Let ttj still denote the smooth global branches of characteristics of A lifted to X. 
Let U{t) be the fundamental solution to p.ll) with operator P lifted to X. Let 
U {t, X, y) be the integral kernel of U(t), x,y E X. Let 

(5.1) U {t, x,y) = ^U (t, x,y), x,y e X, 

py=y 

where x is any point of X such that px = x and the summation is carried out over 
all preimages y of y. Since (15. Ij) is invariant under permutations of the leaves of the 
cover X, the kernel U {t, x, y) is independent of the choice of x. Equation and Cauchy 
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data are automatically satisfied, so ()5.1|) gives a fundamental solution to the system 
(HH) on X. 

The following Lemma gives some estimates for the operator norm from to for 
Fourier integral operators in terms of L°° norms of the amplitude and its derivatives. 
Recall that by we denote the operator norm of T from to H^. 

Lemma 5.2. Let T he a Fourier integral operator 

(5.2) Tu{x)= [ [ e^^(^'^'«)a(x,2/,0«(2/)«, 

where X is an open set in MJ^. Assume that a G S'^ , q & is an amplitude of order q 
and has compact support with respect to x and y in X . Assume also that dip/dx 7^ 
for 7^ 0. Then T extends to a hounded operator from L'^{X) to H^'^^'^^^{X) with 

(5.3) ||r||_,_„_i <C||a(0^'|| Q<\q+n + l\ , 

where {0 = (1 + leH'^'- 

Moreover, let us assume in addition that the support of a helongs to a conical set 
with respect to ^ which does not depend on x, i.e. 

supp a C Hi = {(x, : a; G X, (y, — ) G S, ^ 7^ 0}, 

where S is suhset of S*X. Then for any 6 > we have 

(5.4) \\T\\-g-3n/2-i-s < C meas(E)\\a{^y^\\c\q+n+i\ , 

where meas is the canonical induced measure on S*X and constant C may depend on 
the size of the support of a with respect to x and y. 

Proof. First we consider the case g + n + 1 < 0. Differentiating ()5.2|1 — (g + n + 1) 
times with respect to x we obtain the integral with amplitude of order — (n + 1). The 
integral with respect to ^ converges absolutely, so T extends as a bounded operator 
from L'^{X) to (7^*^'^+"+^) with estimate ()5.3|) . Estimate ()5.4|) clearly follows from this 
as well. 

We will now consider the case g + n + 1 > 0. For a smooth function f , let us 
consider the bilinear form {Tu,v). Let us define operator L as the transpose of *L = 
(1 + \dxip\^)^^{l — id^ip ■ dx). Since dip/dx 7^ 0, integrating by parts g + n + 1 times 
with operator L, we obtain an absolutely convergent integral with respect to ^, and 
the estimate 

\{TU,V)\ < C\\u\\L2\\v\\q+n+l\\a{C)''^\\ci+^+^- 

This implies (j5.3|) . 

We can slightly modify this argument to obtain an estimate of operator T acting 
from L°°(X). Indeed, 



\{Tu,v)\ < 



e 



i^(x,j/,o^g+n+i ^^(^^ i)v{x)) u{y)dydidx 



where 



/ I^{x)dx, Ii{x) = / \a{x,y,C)\d^dy, 

Jx Jx JR" 
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with some amplitude a of order —n — 1. We can estimate a by {^) " and from the 
embedding theorems it follows that 

\(Tu,v)\ < C\\v\\q+n+i\\u\\n/2+s^eas(E)\\a{0~'^\\c<i+r^+i, 
implying estimate ()5.4j) for the norm of operator T acting on L^(X). □ 

The following theorem shows that if a smooth function on a bounded interval has 
zeros only of finite order, then the measure of the set where this function is small is 
also small. Moreover, if we have a family of such functions continuously dependent 
on a parameter, we can estimate measures of sets where functions are small uniformly 
for all parameters varying over compact sets. 

Theorem 5.3. Let W (Z MJ^ be compact and let < T < oo. Let a real valued 
function f = f(t,w) be continuous in w E W and smooth in t E [0,T] up to the 
boundary of [0,T]. Let M E N and suppose that for each w E W function f{-,w) 
has zeros with respect to t of order not greater than M. Then there exist C > and 
Co > such that for all e < tQ we have 

sup meas{t G [0,T] : \f{t,w)\ < e} < Ct^'"^^'' . 

Proof. For e > let 

S(w;,e) = {tG[0,T]: \f{t,w)\<Ce}. 

Let K{w) be the number of zeros of function f{-,w) with respect to t G [0, T] and 
let K be the maximum of K{w) over w G W. It is obvious that K is a finite number 
due to the condition on zeroes of / and compactness of W. 

Let Qfp > 0, p G N, be a decreasing sequence of positive numbers which we will 
choose later. Let us define sets SP(w,e) by setting 

(5.5) («;, e) = {t G [0, T] : |/(t, w)\<Ce,..., 

\dr'fit,w)\ < Ce^^-\ \dU{t,w)\ > Ce^^}, peN. 

We claim now that there exists eo > such that for all < e < eo, and all w E W, 
we have e) = for all p > M. Indeed, if it is not so, then due to compactness 

of W there are sequences t„, Wn and e„, converging to t* G [0,T], w* E W and zero, 
respectively, such that 

and consequently 

dUit*,w*) = 0, p = 0,...,M + l, 

which is impossible. It follows now that the set S(w,e) may be presented as the 
following union of sets 

M 

(5.6) E{w,e) = [jEP{w,e). 

p=i 

The idea of the proof now is to show first that the number of connected components 
of sets Tj^{w, e) is finite and can be estimated uniformly over all w and e. Then we 
will show that the size of each connected component is small and can be estimated by 
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gi/2M^ which will imply Theorem 15.31 These statements are proved in the following 
two lemma. 

Lemma 5.4. There exists eg > such that for all < e < tQ the inequality 

(5.7) A{EP{w, e)) < K{M + 1)^, Vw G 1 < J9 < M, 
holds, where A(S^) is the number of connected components of TP. 

Proof. For simplicity let us first consider the case p = 1. Let us assume that ()5.7|) 
is not valid. Then there exist sequences Wn converging to some w* & W and 
converging to such that 

(5.8) A{T}{w^,tn))>K{M + lf. 

Let us now choose small enough ei > such that each connected interval in the closure 
of S(w*,ei) will include one and only one zero of function f{t,w*) with respect to t 
and will not include zeros of derivative of dtf{t,w*) different from zeros of f{t,w*). 
This is possible because if we had an infinite number of zeros of dtf{t, w*) approaching 
a zero t* of f{t,w*), it would mean that dtf{t*,w*) = and that in fact t* is zero 
of d^f{t,w*) for all A; > 1, which is impossible since we assumed that all zeros of 
f{t,w*) are of finite order not exceeding M. It follows now that A(E(w*,ei)) < K. 
Since w„ converges to w*, by continuity we also have 

(5.9) E\wn,en) cT{w*,ei), 

for sifficiently large n. From (j5.8j) and (j5.9p it follows that there exist (M + 1)^ + 
1 > 2M + 1 connected components of S^(tu„,e„) which are all contained in one 
of the connected components of closure of T,{w*,ei). Let us denote this connected 
component of S(w*, ei) by /. 

According to the definition of sets S^(wji,e„), between two of its connected com- 
ponents function f(t,Wn) must become relatively large (i.e. > Ce) or its derivative 
dtf{t,Wn) must become relatively small (i.e. < Ce°'^). From this we see that in 
the first case dtf{-,Wn) must become zero at some point between these components, 
while in the second case d^f{-,Wn) must become zero at some point there. Since the 
number of components of E^(w„, e„) in / is at least 2M + 1, it follows that at least 
one of these two cases occurs at least M + 1 times. 

Let us consider these cases separately. In the first case, for sufficiently large n, 
we have at least M zeros s;'^ < . . . < of the derivative dtf{t, Wn) contained in /. 
It follows that dff{-,Wn) has at least M — 1 different zeros in J, dff{-,Wn) has at 
least M — 2 different zeros in /, etc. In particular, there are G [s^, s^] such that 
af/(r„^w„) = Oforfc = l,...,M. 

Using compactness of I and continuity of /, it follows that there are subsequences 
s\_, . . . , s^'_, and Wn converging to s];, . . . , s^^ and r^, respectively, which are all 
contained in J, and Wn — > w* . Since functions f{t,w) are smooth in t we have 

dtfisl w*) = 0,..., dtfis!f, w*) = 0, and 9f /(r^ w*) = 0, = 1, .., M. 

Since there are no zeros of derivative dtf{-,w*) on the interval I except may be a 
point t* which is zero of function f{-,w*), we obtain si = . . . = s^^ = t*. From 
G [s^, s^^] it follows that = t* for all = 1, . . . , M, which means 

d^f{t\ w,) = for alU = 0, . . . , M. 
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But this is impossible since function /(-, w*) may have zeros only of order M in view 
of our assumption. 

The second case is sightly different. Here, for sufficiently large n, we have at 
least M zeros < . . . < of the second order derivative dff(t,Wn) contained in 
/\E^(w„,e„). Because in this case we assumed that S^(w„,en) breaks into at least 
M + 1 components due to the failure of condition \dtf{t,w)\ > Ce°% it follows that 



As n tends to infinity, we can choose subsequences of converging to some G /. 
Because dtf{-,w*) does not have zeros in I except may be for some t* E I which is 
also the unique zero of /(■, w*) in J, we get that si = . . . = s^^ = t*. From ()5.10|) we 
also have dtf(t*,w*) = 0. 

Now, to deal with higher order derivatives of / at t*, similar to the first case, we get 
a collection of G [s.^, sff] such that 9f /(r^, Wn) = for A; = 2, . . . , M + 1. Because 
of compactness has a subsequence converging to some , k = 2, . . . , M + 1. Again, 
we must have = t* and hence also d^f{t*, w*) = 0, for all A; = 2, . . . , M + 1. Since 
we already showed that f{t*,w*) = dtf(t*,w*) = 0, this contradicts the assumption 
that f{-,w*) may have zeros of order up to M. 

The argument for p > 2 is similar. For t between two connected components of 
E^, at least one of conditions in ()5.5|1 breaks down. Note that since we assumed 
that the total number of components of e^) is larger than K{M + 1)^ and the 

number of components of the larger set ei) is at most K, we will have at least 

(M + 1)^ + 1 components of S^(i/7„, e„) in I. This means that these p + 1 conditions 
fail at least (M + 1)^ times. Since p + 1 < M + 1, there is a condition that will fail 
at least M + 1 times. 

In the case this is the last condition \dff{t,w)\ > Ce"*" that fails while conditions 
\dlf{t,w)\ < Ce"* remain valid, we can argue similar to the second case of the 
argument with p = 1. In this we need to have at least M — p different zeros of 
df~^^ f{t,w), which is the case if we have at least M — p + 1 such components in /. 
This is true since we have at least M + 1 such components. 

If one of the other conditions fails, let us take the smallest i for which condition 
\dlf{t, w)\ < Ce"* fails M + 1 times. Again, we need to have at least M — i different 
zeros of dl'^^ f{t,w), which would follow from having M — z + 1 such components in 
/. This is again true since we have at least M + 1 such components. □ 

Our next step is to show that each connected component of S^(w,e), for p = 
1, . . . , M, is small enough. 

Lemma 5.5. The length of each connected component of T,^{w, e), for p = 1, . . . , M, 

is no greater than Ce°'''~^~"^ . 

Proof. Let / be a connected component of T,p{w, e) and let t* G /. We are going to 
estimate the maximal shift Sq > such that t* + 6 E I for all < 6 < 6o- Recalling 
the definition of T,p{w, e), we see that 



(5.10) 



dJ{slw,,)\<Ce"^\ I 



1,...,M. 



(5.11) 



\dr'f{t* + 5,w)\<Ce' 



■a. 



■p-l 



|af/(r + 5,u;)| >Ce"^ 
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where t** is some point between t* and t* + 6. Since t*,t* + 6 G /, it follows that 
t** e I and 

\dff{t*\w)\\S\ < \dr'fit* + 5,w)\ + \dr'f{t\w)\<2Ce'^^-\ 

From this and ()5.1H) we obtain that \6o\ < Ce°'p~^~°'p. Consequently, the length of 
each connected component of SP(w,e) can be estimated by Ce°'p^^~°''' . □ 

Now we can finish the proof of Theorem 15.31 Let us choose ak = 1 — k/2M, 
k = 0, . . . , M. According to Lemma 15.51 the length of each connected component 
of TP{w,e) can be estimated by Ce^^^^. Then according to Lemma l5.4| the size of 
SP(w,e) can be estimated by CK{M + l)^e^/^^. Because of decomposition ()5.6p the 
size of S(w, e) is estimated by Ce^''^^'^. Statement of Theorem l5.3l is now a consequence 
of continuity of / with respect to w and compactness of W. □ 

The following interpolation lemma shows that if a function u can be decomposed 
for all sufficiently small e into a sum of u^i'' + with a good estimate for the norm 
of in a "bad" Sobolev space with small index p, and with a bad estimate for 
the norm of in a "good" Sobolev space with large index r, then u belongs to 
some intermediate space H'^ with p < q < r. 

Lemma 5.6. Let u G H^{W^). Suppose that for all small enough e there is a repre- 
sentation u = u^i^ + U2'' such that 

(5.12) Wu'f^ : H^'W <Ce-^, Wu^^'' ■ H^'W < Ce^ , p < r; S,T > 0. 
Then u G if«(M") jor any 

(5.13) q < {pT + rS){T + S)-\ 

Proof. Let a sequence {aj} be such that = ao < cti < . . . and aj +00 as j —>■ +00, 
and let a sequence {bj} of positive numbers bj > tend to zero. Let us assume that 
u G H'^{W^) for some q. Then 



\u 



+00 / „ , — ^ „ _ — 



j=0 \"'aj<l5l<«j+l >^aj<l?l<aj+l 
+00 / „ ^ — _^ 

< 2^ max{(a,)^(«-^), {a.^^f''-^^} / \uf^^\\l + m'd^ 

+ max{(a,)^(^-'-), {a,^^f'^~^^} [ l^'^Wl + m^d^] . 



It follows from this estimate and (j5.12p that 
(5.14) 



\u ■.H^f<CY, (max{(a,)^(^-^), (a,+i)^(''-^)}6f + max{(a,)'(^-''\ {a^^,f^'^-'^]b^^^] 



j=0 



24 



Now we are going to demonstrate that under hypothesis ()5.13|) we may choose se- 
quences {ttj} and {bj} in such way that the right hand side of (j5.14j) will be finite, 
so that conclusion of Lemma will follow. Let us set aj = and bj = with 
some a, (3 > 0. Then the series in ()5.14|) will converge if the following inequalities are 
fulfilled 

2a{q-p) -2PS < -I, 2a(g - r) + 2/3T < -1. 
These inequalities can be transformed into 

a(2T(g - p) + 2S{q - r)) < -T - S, 2(3N > a2{q - p) + I, 
which hold with positive constants a and f3 if and only if ()5.13|) is valid. □ 

6. Spectral asymptotics 

In this section we will prove Theorem 12.51 Let X be a smooth compact mani- 
fold without boundary of dimension > 3. Then an elliptic operator P{x, D) has 
a collection of eigenf unctions and eigenvalues Xj oo. We will be interested in 
distribution of eigenvalues and will find the asymptotics for the spectral function 
A^(A) = tl{j : Xj < A}. One of the most effective methods to study such asymptotics 
is to use an explicit representation for the fundamental solution of the corresponding 
hyperbolic problem We will follow the method developed in [Z|, |2S|, etc., 

for scalar operators, and in [21] for the case M = 1 in Condition C. The following 
proposition was implicitely proved in and formulated in [21j . 

Proposition 6.1. Let X\-,X2 ^ C'^(]R) he such that Xi(0) = 1 and ^ suppx2- 

Suppose that 

Tr JTi^(X2(t)f/(t)) = o(m"-1), ^ oo. 

Then N{X) = ci^-^" + C2{n - ly^X""-^ + o(A"-i). 

In previous sections, and in particular in the proof of Theorem 12.21 we have repre- 
sented the propagator U {t) in (jl.ip as an infinite sum of some extensions of Fourier 
integral operators, i.e. 

oo 

(6.2) U{t) = Y,e-'^'Q„ 

where Qq = I and Tr Qi = 0. So we also have 

(6.3) Tre-^^*gi = 0. 

Let us first consider the contribution of the first term, which is e"*"^*. It is the propa- 
gator for a block-diagonal system, so the asymptotic behaviour of JF^l!^(xo-(^)e~*^*), 
a = 1,2, determines the spectral distribution for a system of independent scalar 
equations. Thus, this is the sum of spectral distributions for scalar operators, which 
are well-known (e.g. [S], [7j). So 

.Fri^Tr(xi(t)e-^*) = Ci/i-^ + c^/.-^ + o(/."-^). 
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We will show that under suitable conditions the asymptotics in ()6.H) are determined 
only by the first term of the series Let us first observe that in view of The- 

orem 12.11 operator J^JLk ^~^^^Qj is compact for sufficiently large K, so it does not 
contribute to the singularity of TrXo-f/, cr = 1,2. Therefore, we can replace U{t) in 
(jfj.lll by a sum of finitely many terms in ()(i.2j) . Also, in view of the second term 

of (|6.2|) does not contribute to (jb.lj) . Other terms of the sum (j6.2j) are of the form 

e-*^* / [ ■ ■ ■[' Z{h)Z{t2) ■ ■ ■ Z{ti+i)dti+i . . . dh, I > 2, 
Jo Jo Jo 

where Z(t) = —ie^^'^Be'^^'^. In this way, after a change of variables si = t — ti, 
S2 = ti — t2, . . . , = ti, we obtain a sum of terms of the form 

IM = Tr y xAsi + ■■■ + s,+i)e'(^i+-+^'+i)^5(s)L(s)ds, a = 1,2. 

where B{s) G L{s) = e-'^n^^^e-'^J^"^ . . . e~'^''+^''+\ and Aj is the j-th block of 
A corresponding to aj. As before, L{s) may be represented as a locally finite sum of 
oscillatory integrals with phase (p{s,x,y,C,), for which 

where Sji,{s,x,p) = ajj,($^^~J ■ ■ Therefore, 

= f [ f xAsi + ■■■ + s,+i)e*(^i+-+^'+i)A'+^^(«.-'.-'.«)6(s, X, X, i)didxds. 
J J X Jr" 

Substituting ^ = firu, r > 0, Icul = 1, we get 

(6.4) /^(/i) = /i" /" [ [ [ Xa{si + --- + si+i)e^'^^''+-+''+^+^^^bT''-^dTdujdxds. 
J Jx is"-i Jo 

To finish the argument for this operator, we can follow |2I] to show that smoothing of 
Q/'s implies ()6.1|) . If we change variables again by p = ^ Sj, sj = puj and introduce 
K = {fij > 0, ^ Kj = 1}, we get 

/ / Xiip)e"'^''^"'^^''"'"'"''"^P^r''~^bdTdpdujdxdK. 

."-1 J~oo Jo 

For fixed u, x, n the point r = —(X] '^kSj,.{pK, x, dip / dx))"^ , p = 1, is a nondegenerate 
stationary point of this oscillatory integral, implying = 0(p"~^^') = o(p"~^). 

This already gives Hormander's first term of A^(A). Now we will show that singu- 
larities at t 7^ do not give essential contributions to the second term of spectral 
asymptotics. In the analysis of Trx2f^, the contribution of the first term of ()6.2|) 
was established in [7]. Let us look at the integral ()6.4|) with respect to si,S2,t with 
fixed X, cj, S3, . . . , S/+1. The stationary phase method with repsect to r, si gives a 
stationary point ip{s,x,x,uj) = 0, r = —{dip/dsi)~^. It is non-degenerate because 
det drdsj^i^Tip) = —ai{s^x^dx^Y ■ This gives the estimate /2(p) = Oi^jji^'^). Since 
diplds2 = only at r = —{dip/ds2)^^ = a2{s,x,dx'^)~^, the phase is stationary with 
respect to S2 only at (x, ^) for which ai{s,x,dx(p) = a2{s,x,dx'^). This is the set of 
measure zero and, therefore, hip) = o(p"~^), which shows ()6.1|) . Finally we note 
that if the support of xi is sufficiently small (i.e. when s is small), terms in ()5.ip . 
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corresponding to x and y in different leaves of X will produce only smoothing oper- 
ators (j5.1|) in view of the finite propagation speed. So their contribution to I{fi) is 
rapidly decreasing and we obtain the result also on X. 
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